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A path integral formulation of the Bianchi I models containing a massless scalar field in loop quan- 
tum cosmology is constructed. Following the strategy used in the homogenous and isotropic case, 
the calculation is extended to the simplest non-isotropic models according to the p and p! scheme. 
It is proved from the path integral angle that the quantum dynamic lacks the full invariance with 
respect to fiducial cell scaling in the p scheme, but it does not in the p! scheme. The investigation 
affirms the equivalence of the canonical approach and the path integral approach in loop quantum 
cosmology. 
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■ I. INTRODUCTION 

Loop quantum gravity [l[ and the spinfoam formalism 0,01] can be considered respectively as canonical quantization 
and covariant quantization of gravity just like what we did in ordinary quantum field theory. Applying the the 
techniques used in LQG to simple cosmology models, we get loop quantum cosmology 4], which is a canonical version 
of quantum cosmology. Naturally, we desire to understand LQC conclusions from a path integral, i.e., covariant 
angle. Starting from the Hilbert space of LQC, there are two kinds of path integrals. One is integral over paths in 
configuration space which leads to a spinfoam like 'vertex summation' rather than a standard path integral 

format. The other is integral over phase space paths 0-[Lll from which we can get a typical continuous path integral 
_r that features continuous paths with weights given by the exponential of the phases. 

In this paper, we will use phase space paths following the methods in the original work of A. Ashtekar, M. Campiglia 
and A. Henderson Q, and extend the LQC path integral formulation from homogenous and isotropic FRW model to 
the simplest homogenous Bianchi type I models [12j . Up to now there are two kinds of 'loop regularization' of the 
gravitational part of the scalar constraint for the Bianchi I models in the literature, referred to as p and p! scheme 
Each of them has its own advantages and drawbacks, pi scheme gives a difference equation in terms of affine 
, variables and therefore the well-developed framework of the spatially flat-isotropic LQC can be straightforwardly 
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adopted. While the p! scheme has better scaling properties 14J, but the difference equation of this scheme is very 
complex. Here we use both p and p,' schemes presented in [l3| to perform a path integral formulation, and to explore 
fSj ' their similarities and differences. 

The organization of this paper is as follows. In the Sec. HH we present the quantum theory of Bianchi I modlcs in 
the volume and connection representations. In Sec. IIII1 we construct the path integral formalism of quantum Bianchi 
I models in phase space using pi scheme, and the construction in p! scheme will be presented in the Sec. IIV1 Finally, 
we discuss the results obtained from those two schemes in the Sec. IVl 



II. LQC OF THE DIAGONAL BIANCHI I MODELS 



Bianchi cosmologies are homogeneous cosmological models, in which there is a foliation of spacetime M = £ x K 
such that £ is space- like and there is a transitive isometry group freely acting on S. Thanks to these symmetry 
properties, in classical theory, after the diagonalizing and rescaling processes of the Ashtekar variables (see jT^] ) , the 
only non- vanishing constraint (scalar constraint) for the Bianchi I cosmology (coupled with a massless scalar field <f) 
is given by 

C = - g7rG , a (c 1 piC 2 p 2 + cVlC 3 p 3 + C 2 p 2 C 3 p 3 ) + (1) 
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where p^ is the momentum of the scalar field, while c 1 and pj are the reduced and rescaled Ashtekar variables. The 
Poisson brackets of c 4 and Pj are 

{c\p J } = ^ 1 G8). (2) 

In the isotropic models, there is a reflection symmetry H(p) = — p, which corresponds to the orientation reversal of 
triads. These are gouge transformations which leave the Hamiltonian constraint invariant. In Bianchi I case, we have 
three reflections IVs corresponding respectively to the reversal of one of the triad while leaving the other two fixed. 
The Hamiltonian flow is invariant under the action of each IT fl4j . Hence it suffices to restrict our attention to one 
of the octant of three p^s. In this paper, we focus on the positive octant in which all three p^s are positive. 

A. p, scheme 

In the canonical quantization theory, i.e., loop quantum cosmology, there exists a comprehensive construction of 
the operator corresponding to the Eq. (TT|) in the full Loop quantum gravity scheme The full LQG quantization 
proposes two kinds of quantum corrections to the scalar constraint, one of which is the inverse triad correction by 
using a so-called Thiemann trick, while the other one is holonomy correction which means to correct the curvature 
term by usin g a SU(2) holonomy. In this paper, we start out setting the lapse function N — IP1M2P3I m the classical 
theory as in [14| . This treatment doesn't only make the Hamiltonian constraint simpler, but also avoids the use of 
inverse triad correction. Hence, we can use the regularized constraint [l3j to replace Eq. (TTJ), in which only the LQG 
effects from holonomy correction are taken into account: 



c = 1 
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sinf/iic 1 ) sin ((J.2C 2 ) sinf/iic 1 ) sin (/13c 3 ) sin (/i2C 2 ) sin (/13c 3 ) 

= -Pi = -P2 H = -Pi = -P3 H = -P2 = L P3 

Mi Mi M2 



f=0. (3) 



In n scheme, /ij = y/ A/\pi\, where A = 4\/ r 3Tr r ylp l is the 'area gap' [Ta|. Usually, we can make an algebraic 
simplification by introducing new phase space variables: 

rA = — T=, (4) 



6ttG l\/p 

where the index i does not sum over and their Poisson bracket is given by 

{&>,} = 2<5j. (5) 
Then the kinematical Hilbert space of the gravitational part is given by 

H-Kin — L 2 (M.Bohr,dHBohr) 3 , (6) 

with the orthonormal basis \vi, ^2, ^3) = Wi) W2) l^) satisfying 

(Vl,V2,V3 \Vl,V2,h) = HuuViSvzfaSvsfo. (7) 

Any state ^(^i, V21 ^3)) G %Kin can be decomposed in the orthonormal basis as 

1*2,1/3)) = X] *(^1^2, ^3) 1^1,^2,^3) • (8) 

To define the constraint operator, we have to define two kinds of operators first. One is the volume operator defined 

as 

V\ui,U2,v 3 ) = 67tG7v / AVi^3| 1/3 1^1,^2,^3) , (9) 

where V = \ViV 2 V 3 \ 1/3 , and V t \ Vi,V2,v 3 ) = GTTGjVAfi \vi,U2,ua), which is the operator corresponding to pi. The 
other is the unitary shift operator 

frfa) = exp (ifk?) |n) = \vi + 2e h) , (10) 
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where £q is related to the 'area gap' by l\ = 7 2 A, where 7 is the Barbero-Immirzi parameter. Replace the function 
sin fiiC 1 in Eq. ([3]) by the operator sin \XiC % using the unitary shift operator defined above, we get 



sm(/j,ic l ) \ui) = ~~ (exp^c 1 ) - exp(-i^c 1 )) 



2i 

Then we fix the factor ordering ambiguity of Eq. (|3|), by defining the operator 

a , , 1 / sinZ2,c l ~ - sin/2, c l \ , , 
Gih) :=- —g^Vi + Vi 1 1 '- 



1 {\ Vi + 2l h)-\vi-2i Q h)). (ff) 



A VA J 

= - 3z^G 7 [{ Vi + i h) \ Vi + 2i K) - 0* - £ h) \ui - 2£ h)} . (12) 

Notice that [0j, 0j]|i/i, v 2 , v 3 ) = for i =^ j. The gravitational part of the scalar constraint operator can be written 
in terms of these operator as 

d9raV "~8^V (6162 + 6263 + 6163) • (13) 
The action of this operator on arbitrary basis element of the kinematical Hilbert space is [l6[ 
C^\ Vl ,v 2 ,v 3 ) 

= %~ [(fi + 1)(^2 + 1) Wi + 2,v 2 + 2, v 3 ) - (y x - \){v 2 + 1) \vi -2,v 2 + 2, v 3 ) 

o 

-fa + 1)(^> - !) I^i + 2, ^2 - 2, u 3 ) + fa - \){v 2 - 1) Wi -2,v 2 - 2, v 3 ) 
+(j/ 2 + l)(i/ 3 + l)\v 1 ,v 2 + 2,v 3 + 2) - (y 2 - 1)^3 + 1)1^2 -2,^ + 2) 
-{v 2 + l)(f 3 - 1) \vi,V2 + 2,^3 - 2) + {v 2 - \){v 3 - 1) \v u v 2 -2,v 3 - 2) 
+(1/1 + l)(z/ 3 + 1) |z/i + 2, V2, v z + 2) - (wi - + 1) \v\ - 2, f 2 , ^3 + 2) 
-(j/i + l)(i/ 3 - 1) \v x + 2, v 2 , v z - 2) + (^ - l)(i/ 3 - 1) \vx - 2, iaj, i/ 3 - 2)] , 

where the l$h terms are omitted for short. From this equation we can see that in the fi scheme, the space of solutions 
to the quantum constraint is very similar to the isotropic case. v\, v 2 , v 3 are discrete variables and they are supported 
on a specific superselections respectively. They can take the values of z/j = (e^ + 2n^)l^h^ where the parameter 
£i G [0, 1], and rii £ N. As mentioned previously, we focus on the positive octant here. This precondition simplify the 
calculation but do not affect the physical meaning. If we consider the general situation, a factor sgn(pi) will appear 
in the expression of 6^. And as a function on phase space, it does not commute with sin /i^c 1 , hence their product 
as operators is not symmetric. It is necessary to realize this point, although our simplification makes its role less 
important. For the details about this situation, see [l7j . 



(14) 



B. // scheme 

We still use the regularized constraint Eq. ([3]) in p! scheme. Only fa changes to fa] in this situation, and Jx[ takes 
the form 



u'- NA a'- |P2|A a'- WA (15) 
Owning to this condition, it's inconvenient to still use i/i as configuration variables. By introducing new variables A^, 

Xi = (4ttG)V3' ( 16 ) 

the algebra of the p! scheme could be much more simplified 14]. The variable conjugate to is 

_ 2^ 

7 (47rG) 2 / 3 ' 1 ' 
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We can prove that the resulting Poisson bracket between k l and Xj is 



{k\X 3 } = 25] 



(18) 



Then the new orthonormal basis is |Ai, A2, A3) in Hf^, and the fundamental volume operator and unitary shift are 
given by 



and 



V|Ai,A 2 ,A 3 ) =2^G 7 VAw|Ai,A 2 ,A 3 ) 

£ h 



E,X 



exp (i^c 1 ) \ Xi) 



\ + 



XjXk 



(19) 



(20) 



Here, v = 2A1A2A3 is the configuration variable related to the volume of the elementary cell, and the factor 2 ensures 
that this v reduces to the v used in the isotropic case. We can also use Ai, Xj, and v as the basic configuration variables 
to substitute for Ai, A2, and A3. The indices i, j, k = 1,2,3 appeared in Eq. p^|) are required to be different with each 
other, and, as usual, the index % does not sum over. Having these operators in hands, we could get the gravitational 
part of the scalar constraint operator. Its action on arbitrary basis element of the new kinematical Hilbert space is 

m 



C 3rav \X u X 2 ,v) 



~2~~ 



(v + 2)V^7+4|Ax, X 2 ,v)2 - (v + 2)^\X u X 2 ,v)+ 



- (v - 2)Vw| Ai , A 2 , »)o + (v - 2) V|u-4|| A x , A 2 , v) 4 



(21) 



where 



|Ai, X 2l v} 4 



v±A v±2 

Ai, A 2 ,w±4 



v± 2 
v±2 



v±4 



1; 



v±4 

•A 1; — — -A 2 ,w±4) + 
v ± 2 



v±2 
v±2 



v 



Ai,A 2 ,w±4 
Ai,A 2 ,w±4 



x v±2 ^ 

Ai, A 2 ,w±4 



Ai,— — - • A 2 ,w±4 
v ±2 



(22) 



and 



|Ai,A 2 ,w) 



v±2 



•Ai, 



v±2 



A 2 ,w 



v±2 



v±2 



Ai 



, v±2 , 
• Ai, X 2 ,v ) + 



v±2 



■X 2 ,v) + 



v±2 



x v±2 x 
Ai, X 2 ,v 



■ Xi,X 2 ,v 
■ X x ,X 2 ,v 



(23) 



As what we did in Eq. (jl"4")l . the £qH terms are omitted here. From Eq. (|2"2"|) and Eq. (|23| . we can see that the 
variable v is fine, and the wave function only involves (v — <i£oh) , v, (v + Moh) terms. It is exactly the same as in the 
isotropic case. On the contrary, the situation for Ai, A2 is much different. We see that they depend only on the value 
of v. This dependence is through fractional factors whose denominator is two or four units bigger or smaller than the 
numerator. This lead to the consequence that the iterative action of the constraint operator derives only to states 
whose quantum numbers Ai are of the form Ai = Lo f A', here A^ is the initial value of Ai and e is a constant number 
that v = e + 4n£ h, n E N. oj e belonging to the set 



a 



n 

m,n£N 



e + 2m 
e + 2n 



(24) 



where fc™ £ N, and z 6 Z if e > 2, while z = when e < 2. The discrete set J7 e is countably infinite and turns out 
to be dense in the positive real line. The proof of this statement see Appendix D3 of [llj]. So, On the one hand, A; 
is superselected in separable sectors, have countable numbers of values. On the other hand, the dependence on the 
fractional factors makes it can take values in entire positive real line. 
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III. PHASE SPACE PATH INTEGRAL IN p SCHEME 



In this section, we calculate the extraction amplitude A{vf,<f>f,Ui,<j)i) for Bianchi I cosmology (here V is short for 
{y\, V21 v sY)- As we know, in the 'timeless' framework of LQC, the whole information of the quantum dynamics is 
encoded in A(pf,<j)f;vi,<f>i), which is the transition amplitude for our phase space path integral just as in ordinary 
quantum mechanics. We can express it as [f| 



A (vf,<t>f,Vi,(()i) = / da (v f 



e " 



of 



(25) 



where the total constraint operator C is composed of two parts: the gravitational part and the scalar field part 

C = C grav + C matt . (26) 

The two operators C 9rav and C matt are commutative, so they can act on their own Hilbert space, 'H^in anc ^ ^Kini 
respectively. Then following the spirit of @, decompose the fictitious evolution into N evolutions of length e = 1/N, 
and insert complete basis in between each factor, we get 



{Vf,<j>f 



ei aC 



UN-I—Vl 



1 {VN, <PN 



e» 



taC 



VN-l,<pN-l) ■ ■ ■ \yii9l 



(27) 



where (v N ,<p N \ = (vf,<f>f \ and \v ,(f>o) = Wi,4>i)- 

Firstly, we need to calculate the n-th term appearing in above expression. Using the Eq. (1261) . it can be stated as 



e'< 



eaC mat 



K) (^1+1,^+1,^+1 



eaC gra 



I/ 1 I/ 2 7/ 3 \ 



(28) 



As usual, the matter Hilbert space is the standard one, "H^in = L 2 (M,, d<j>). By an ordinary quantum mechanics like 
calculation, the scalar field factor is 



en 



EOC" 











n) = 









2tt 



(29) 



and the gravitational factor is 

/ 1 2 3 

\^n+l, v n+li U n+1 



en 



5 < + ^J< + l.< 5 <+l,< + i £a (^+l^n+l^n+l C 9 ™" V li V li V l) +°0 2 ), 



(30) 



where we made an expansion in e, and C(e 2 ) is the high order terms of the expansion. Using Eq. (J3|), the matrix 
element of C grav is 



/ 1 2 3 

\ v n+li ^n+li v n+l 

9ttG 



C gr 



I 9 3 \ 



v n+l + v n v n+l + v n , x x x x x x 



. ^n+l + V l V l+1 + / c c c c 

H 2 2 V "n+l> 1 ^+ 2 "n+l> 1 '»+ 2 "i+l^n _ " " ' 



-2 c, i/J + 1 ,^+2 u ^ j 1 i + 1 ,i/i 



"^+1.^+2^ + 1.^-2^ + !.^ + ^+i»»'5-2^ +1 ,i/8-2^ +1 ,j/1 



^n+1 + ^n+l + 



(31) 
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As demonstrated before, we choose a specific superselection of vi such that z/j = 2rii£oh, where n, € N. Then we can 
use the Fourier expansion formula of the Kronecker delta 



x i tu 
ov v = — 



dbe- lb ("'-") /2h , 



(32) 



where, obviously, b is the conjugate variable to v and take values in the range (0, tt/Iq). By using Eq. (|32p . then Eq. 
f[3T| can be expressed as 



/ 1 2 3 



i 9 3 \ 



— ) / db n+ ie 



. e r : 



^ ^.4-1 + ^ 



— sin i bl +1 sin £ &„+i 



- sin^o&n+i sin ^o^+i 



*4+i + K v n+1 



-^sin^+isin^ 3 n+ i+0(e 2 ) 



(33) 



Now according to the calculation before, we could put the scalar field part and gravitational part together, and Eq. 
(l27l) takes the form 



(v f ,<j>f \e laC ViAi) 

3N 



dbjy ■ ■ ■ dbi ■ [ — 



:SA r 



dp N ---d Pl e* SN +0(e z ) 



(34) 



where 



5. 



AT 



E/ 0n+l — 0n ^n+l ^n+l — V n &n+l ^n+l — 1/2 ^n+1 v n+l 
[ Pn+1 ■ 



ri=0 



+a 
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Pn ~ 



K+i + K K+i 



K, 



sin^o6^ +1 sin£ b n+l 



K <4-i + v 



'n+1 "n "n+1 



- sin £ob^ +1 sin£ &: 



3 

n+1 



vl ^4-1 + V? 



'n+1 1 "n 



-sin46„ +1 sin^ 6: 



3 

n+1 



(35) 



This is the 'discrete-time action' of the Bianchi I cosmology. Just like in |{|, the final step is to take the limit N — » oo, 
and because the variable i/j is discrete, it is impossible to interpret the {v n+ i — v l n )/e as a derivative, we should 
transform the terms (f„ +1 — f n )/e into {b l n+1 — b l n )/e: 



N-l 



\ (vNbN - Voh) = \ ^2 ( b n+lV n +l ~ V n b n ) . 



n=0 



Then 



(36) 



JV-l 

E 



b n +l Vn+\ — V n 



N-l 

71=0 



Vn b 



n+1 — "n 



+ -^(h^Q - b N v N ). 



Using this transformation, Sn can be rewritten as 

-l-K . Vn b n+l-t> n , vlb 2 n+l -b 2 n v\b\ 



Sjy — 



- b 3 
n+1 u n 



+a 



Pn - 2 nG 



'n+1 T % " n+ i -r " n 



sm£ bl +1 sme bl +1 



v l+l + v n v l+\ + "n . , 2 . „ , 3 . V n +1 + V n v l+l + ^ 

o 7> sm46 n+1 sm46„ +1 + ; 



(37) 



- sm£ bn +1 sm£ b n+1 



-(vfbf - vA). 



(38) 
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Finally, take the limit N oo, we have 

A(p h $ f -n,fa) = / a I [Diy q (T)}[Db q (T)}[Dp(T)}[D0(T)}e*», (39) 



where 



s ' = / <It <J p(j) + —v -b — a 



g 

p 2 - -nG^u 2 sinlob 1 sin^fe 2 + u 2 ^ sin£ 6 2 sin£ & 3 



(40) 



As in the isotropic case, we use the subscript q here to emphasize that now we can take a sum in the geometrical 
sector including only the 'quantum paths'. Also, we use the same trick to transform it to the familiar format 

^ rir/tta poo poo 

— V / db n -> dv n / db n . (41) 

*0 v Jo J -oo J -oo 

Using Eq. (|41[) . we get the final expression of the extraction amplitude 

Apfrffiifuk) = J da J [Du(T)}[Db(r)}[Dp(r)}[D^T)}e^ s , (42) 

where 



S = J dr l^pcj) — —b ■ v — a [p 2 



-^ixGiv^v 2 sin^ob 1 sin£ 6 2 + v 2 v z sin£ 6 2 sin^fo 3 

+^V 3 sin46 1 sin46 3 )]}. (43) 

This is the path integral formulation and its action we desired. Because we used Eq. (|41[) . the integration about Ui 
and bi are taken from — oo to oo, with the boundary terms in Eq. (t40|) canceling each other out. And then we are 
able to integrate over all paths in the classical phase space as in usual path integrals. 



IV. PHASE SPACE PATH INTEGRAL IN p! SCHEME 

The construction of a phase space path integral in p,' scheme is more complex and subtler than in p scheme. As we 
illustrated in III B[ the new orthonormal basis in H 9 ^^ can be |Ai, A2, A3) or |Aj, Aj, v). In order to compare with the 
p scheme situation conveniently, we choose the former one. Following the procedure employed in IIII1 we are going to 

calculate the extraction amplitude A (Xf,<f>f,\i,(j)ij. 

The matter part of the extraction amplitude we use here is the same as in p scheme, so we ignore it for a little 
while and directly use it at the end. Decompose the fictitious evolution into N evolutions of length e = i, the n-th 
term of the gravitational part is 

= <S(Ai l+1 - A")(5(A™ +1 - AJ)5(A 3 1+1 - A3 1 ) 

+ l -ea{Xl +1 1 Xi +1 ,\l +1 \C arav \Xl,\ , i,X^)+0{e 2 ). (44) 

As we known, although the number of possible values of \ is countable, there does not exist an obvious superselection 
with respect to it, therefore, we have to work with the entire positive real line it spans. So, for the possibility to perform 
a path integral, in a sense, we have used an approximation here. The Dirac delta is used rather than Kronecker delta 
to express the inner product <A" +1 |A"), which means we will take \ as continuous variables below. This assumption 
may generate a little confusion. When we think Xi is continuous, the quantum theory will be changed which caused by 
the changing of the Hilbert space. But, what we interest in the following is the effective action, which is a semiclassical 
quantity itself. And consider the properties of Xi, it is reasonable to perform a semiclassical approximation that let 



e,r, 



eaC 1 
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Xi to be continuous in the following calculations and do not affect the quantum theory. It is worth to emphasize that 
we will focus on the semiclassical situation from now on. 

Then, make a algebraic transformation: v 1 = v + 2, = v'Xi/(v' — 2); v" — v — 2, A" = v"\i/(v" + 2). In term of 
Eqs. (JUJ), and the matrix element of C 9rav is expressed as 



\ A 1 ; A 2 ' A 3 



ttG 



(45) 



where 



(A) = ±(A? +1 ,A2 +1 ,A£ +1 



_i_/\n+l \n+l \n+l 



1 \/n 



\/n \n\ _i_ /\n+l \7l+l \n+l 
n> A 2 ) A 3 y ± \ A 1 > A 2 ) A 3 



\n \n _i_ A)^ \'n\ i /\ri+l \n+l \ra+l 
2 A^A^" 3 / \ 1 ' 2 ' 3 



\n_i_ , n wn 

1 \«\/n' A 2' A 3 
A 2 A 3 

12 AfAg 



±(A? +1 ,AJ +1 J AJ +1 



A l i A 2i A 3 ^ 



t h 

A i" A 2 



±<A^+ 1 ,A' 2 l+1 ,A^ +1 



A l ) A 2 ) A 3 



e h 

a?a 2 " 1 ' 



(46) 



and 



±(A? +1 ) A5 +1 ,A^ +1 



A 2 A 3 



A 2 '" ) AA±(A? +1 ,A^ +1 ) A^ +i 



A l T w,„ , A 2 , A 3 
A 2 A 3 



±(A? +1 ,A5 +1 J AJ +1 



12 A^™ 3 



A l ? A 2 j A 3 T 



3 1 \tln \ n 
A l A 2 



Employ the identity 



±(a; i+1 ,a 2 i+1 ,a^ +1 



\n+l \n+l 



\ //n \ 71 -,- 

A l i A 2 + 



\7l 

//7l \7l ' 3 



A'/" A 



±(A? +1 ,A 



2 i A 3 



\ n \ I In \ n 
A l > A 2 i A 3 + 



A?A 2 '" 



(47) 



5 (A' - A) = 



-Hoc 



dfce -«(V-A)/2» 



(48) 



where k is the variable conjugated to A as we stated in Eq. ([T71) . we are led to the following expression for Eq. (j44 



(X" +1 



eaC sr 



A' 



1 



x 



(*/)" ( e l (- fc ? A ^- fc ?^^)^ _ e i(fc?Ar^-fcJAJ^-)«, . . . + cydic terms ^ 

+ (v") m (e< k " x "^+ k ^ x ^) io - e K-k?K^+^y^)t . . . + cycHc terms ^ + G ( e 2) I _ ( 49) 

Let £ k™X"/v n = a iy £ k"Xf/ («')" = a' it and iokfXf/ (v") n — a" for short. The exponential terms of this equation 
can be expressed as 



v {[cos(a4 + 02) — i sm(a[ + 02) — cos(a' 1 — 02) — i sin^ — 02)] + cyclic terms} 
+v" {[cos(a'/ + 02) + i sin(a" + 02) — cos(a'/ — 02) + i sin(a'/ — 02)] + cyclic terms} . 



(50) 



Here, we introduce another approximation: assume the volume of the elementary cell is large enough, then we can let 
v = v' — v", and therefore the i sin terms in Eq. ([50|) are canceled with each other. So we have 



X n+1 
1 



A' 



(7T)3 

x I sin 



7,71 \ 71 h n \ n 

""1 A l n ■ 2 A 2 

1 sm 



<u + sin + sin sin ^ i + G(e 2 ) 



v r, 



(51) 
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Combining the scalar field part presented in Eq. (|29|) . the extraction amplitude takes the form 



A(Xf,<j> f ;XiAi) = (\ f ,4>f 



y — 



2/V 



,, , dk N ■■■dk 1 -[ — ) I dp N ---d Pl e^ N +0(e 2 ), 



where 



N-l 

kfX 



x sin 



e 2 e 

fc£A 

v 



{p n f- 4nGVv n v n + 1 v r 



v" 



i k n X r 
-to sin -2-^-l Q + sin -^-£ sin ^— M + sin :lL;lL 4 sin ^— M Q 



/t 3 A3 



h.n\n 



/t 3 A3 



Finally, taking N — > oo, we get 



DX(t)] \d%(t) \ [Dp(r)] [D4>{t)] e* s 



where 



5 



r/r 



-fc • A - a [p 2 - 4ttG 



x« ( sin «o sin to + sin 



fc 2 A 



feA 2 . fc3A 3 „ 
to sin «o 



. fciAi. . fc3A 3 „ 
■ sin to sin .eg 



(52) 



(53) 



(54) 



(55) 



Eq. (|54"|) and Eq. ([S3)l are the final expressions of the extraction amplitude and its action in // scheme. Notice 
that Xi has already been taken as a continuous variable, therefore, it is reasonable to interpret the ( A™ +1 — A™ ) /e as 



derivative here, and the integration ranges of A and k are taken from — oo to oo directly. We no longer need to play 
the trick which was used in p, scheme and the isotropic case to transform the integration range. Integration is over 
all trajectories in the classical phase space from the beginning. 



V. DISCUSSION 



In this paper we extended the phase space path integral formulation of Friedmann space-times to anisotropic Bianchi 
I models, and performed the calculation both in p, and p! scheme. We restricted the matter source to be a massless 
scalar field and focused on the positive octant in which all three pi are positive. 

Since pi = \j A/ P i in p scheme, by setting ~ p 3 / 2 '', just as v ~ p 3 ^ 2 in the isotropic case, the procedure used in 
@ could be implemented directly. The formulation of the extraction amplitude then resembles three copies of the 
isotropic model as we desired. To incorporate the p! scheme, we have to overcome a few obstacles. The main problem 
comes from the fact that the expression of the p!^ which was given by Eq. (1151) . is much more complicated in p! scheme 
than in p scheme. This leads to the consequence that we have to make an algebraic simplification by introducing 

1/2 

Xi ~ Pi ' , and then the wave function is dragged along the Ai direction by the unitary shift operator E\. However, 
the affine distance involved in this dagging depends on A2 , A3 . From the difference equation of the gravitational part 
of the total constraint, it can be seen that the new variable Ai has a discrete spectrum but does not support on a 
specific Ai-lattice. We can not write an integral representation for the Kronecker delta by using a standard residue 
calculation like Eq. (f3"2"j) . 

In order to perform a path integral formulation in p[ case, we have used two approximations in HVl The first one 
is taking Ai as continuous variable, then one can use the Dirac delta to express the inner product. The second one 
is to presume the volume of elementary cell is large enough such that (v — 2) ps v ps (v + 2). Actually they are the 
same thing. In general, our treatment presupposed that the the calculation is away from the Planck regime, in other 
words, a semiclassical approximation was taken. This is the key point where p' is different from p and the isotropic 
case. In the two later situations the calculation is exact and does not need any additional assumption. 
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If we compare the two path integrals in p and p! scheme, it is obviously that the actions have the same formulation 
as in isotropic case: a classical action plus a quantum gravity correction term, i.e., the sin term. Ignoring the constant 
terms, from Eq. fj43]) and Eq. ([55]) . one can write down directly the two quantum gravity correction terms in the two 
schemes (denoted by T and T' respectively). In p scheme, 

T = ^^sinlo^smiob 2 + is 2 v 3 sm£ b 2 sm£ b 3 

+i/Vsin£ b 1 sin£ 6 3 . (56) 

In p' scheme, 

rr> 2 f- Wi ( . k 2 \ 2 . fc 2 A 2 . fc 3 A 3 
I = v \ sin lo sin to + sm to sm 1$ 

\ V V V V 

. . feiAi . fc 3 A 3 . \ 

+ sm £ sm 1 . (57) 

w v J 

These terms do not look the same. By using the explicit expressions of bi,Xi and hi, T and T" could be expressed as 
T = v v sin/lici sin/i2C2 + ^ v sm P2C2 sm //3C3 + v v sin/iici sin/i 3 c 3 , (58) 

and 

T' = v 2 (sin/2'jCi sin/2' 2 C2 + smp,' 2 C2 sin/2 3 c 3 + sin/i^ci sin/IgCs) . (59) 

Notice that ~ p^ 2 and i> ~ (piP2P3) 1 ^ 2 - Now we find that T and T' are extremely similar in the formulation. The 
difference is the immediate cause for the fiducial cell dependence on the semiclassical limit in p, scheme. Consider the 
fiducial cell V. Its volume is Vb = L1L2L3. Notice that Ci ~ Lihi and pi ~ LjLkCijak, so if we rescale V as 

= L X L 2 L 3 -»■ Vq' = hL x hUhL z = hhhVo, (60) 

then 

/2ici -^= /2iCi, (61) 

V '2'3 

and 

/zici ->• /iici. (62) 

This is the reason why /I' scheme has better scaling properties. These properties are already known in canonical 
quantization program. In path integral formulation, we have proved again that the quantum dynamic does depend 
on the choice of the fiducial cell in p, scheme, while it does not in pi scheme. Furthermore, if we take p\ = p 2 = Pa, 
and therefore Vi = v, then we find T = T 1 ~ v 2 sin 2 pc. This is precisely the quantum correction term in Friedmann 
model. Hence, it is rational to forecast that in the path integral formulation of Bianchi I models the quantum bounce 
will also replace the big bang singularity due to the quantum gravity correction term. All of these properties provide 
a strong evidence for the equivalence of the canonical approach and the path integral approach in LQC. 
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